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Abstract 

We theoretically propose quantum spin pumping mediated by magnons, 
under a time-dependent transverse magnetic field, at the interface between 
a ferromagnetic insulator and a non-magnetic metal. The generation of 
a spin current under a thermal equilibrium condition is discussed by cal- 
culating the spin relaxation torque, which breaks the spin conservation 
law for conduction electrons and operates the coherent magnon state. 
Localized spins lose spin angular momentum by emitting magnons and 
conduction electrons flip from down to up by absorbing the momentum. 
The spin relaxation torque has a resonance structure as a function of the 
angular frequency of the applied transverse field. This fact is useful to 
enhance the spin pumping effect induced by quantum fluctuations. We 
also discuss the distinction between our quantum spin pumping theory 
and the one proposed by Tserkovnyak et al. 



1 Introduction 

A new branch of physics called spintronics has seen a rapid development over 
the last decades. The central theme is the active manipulation of spin degrees 
of freedom as well as charge ones of electrons; spintronics avoids the dissipa- 
tion from Joule heating by replacing charge currents with spin currents. Thus 
estabhshing methods for the generation of a spin current is significant from view- 
points of fundamental science and potential applications to green information 
and communication technologies. [T] 

A standard way to generate a (pure) spin current is the spin pumping effect 
at the interface between a ferromagnetic material and a non-magnetic metal. 
There the precession of the magnetization induces a spin current pumped into 
a non-magnetic metal, which is proportional to the rate of precession of the 
magnetic moment. The precessing ferromagnet acts as a source of spin angular 
momentum; spin battery. This method has been developed by Tserkovnyak et 
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al. [5] Now their formula has been widely used for interpreting vast experimental 
results, despite the phenomenological treatment of spin-flip scattering processes. 

Thus we reformulate the spin pumping theory, through the Schwingcr-Keldysh 
formalism, to explicitly describe the spin-flip processes. A generation of the 
pumped spin current is discussed on the basis of the spin continuity equation 
for conduction electrons. Moreover the utilization of magnons, which are the 
quantized collective motions of localized spins, has recently been attracting 
considerable interest. [31 SI [17] Hence, the focus of the present work lies also 
on the contribution of magnons to spin pumping. We treat localized spins as 
not classical variables but magnon degrees of freedom. Therefore we can cap- 
ture the (non-equilibrium) spin-flip dynamics, where spin angular momentum 
is exchanged between conduction electrons and localized spins, on the basis of 
the rigorous quantum mechanical theory. Consequently, we can reveal the sig- 
nificance of time-dependent transverse magnetic fields, which act as quantum 
fluctuations, for the spin pumping effect. 

We consider a ferromagnetic insulator and non-magnetic metal junction|3] 
shown in Fig. [T] where conduction electrons couple with localized spins S(x,t), 
X = {x, y, z) G R^, at the interface; 

nc^^-2Jal dx S(x,i) • s(x,t). (1) 

xG (interface) 

The exchange coupling constant reads 2 J, and the lattice constant of the fcrro- 
magnet is og. The magnitude of the interaction is supposed to be constant and 
we adopt the continuous limit in the present study. Conduction electron spin 
variables are represented as 

■nX=t.i (2) 
= (cta^c)/2, 

where cr-' are the 2x2 Pauli matrices; [cr-', cr'^] = 2iejki<j'', {j,k,l — x,y,z). 
Operators c^/c arc creation/annihilation operators for conduction electrons, 
which satisfy the (fermionic) anticommutation relation; {c,j(x, t), c^(x', t)} = 

We focus on the dynamics at the interface where spin angular momentum 
is exchanged between conduction electrons and the fcrromagnct. We suppose 
the uniform magnetization and thus localized spin degrees of freedom can be 
mapped into magnon ones via the Holstein-Primakoff transformation; 



5+(x,t) EE 5^(x,t) +iS'^(x,t) 

= \/2§[l - a^{x,t)a{x,t)/{4S)]a{x,t) + 0{S~^^^), 

= V2l'a^(x,t)[l - at (x,t)a(x,i)/(45)] +0(5-3/2), 



(3) 
(4) 
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Figure 1: (Color online) A sehematie pieture of quantum spin pumping mediated 
by magnons. Spheres represent magnons and those with arrows are conduction 
electrons, (a) Localized spins lose spin angular momentum by emitting magnons 
and conduction electrons flip from down to up by absorbing the momentum, 
(b) An enlarged schematic picture of the interface. The interface is defined as 
an effective area where the Fermi gas (conduction electrons) and the Bose gas 
(magnons) coexist to interact; J 7^ 0. Conduction electrons cannot enter the 
ferromagnet, which is an insulator. 



^^(x,^) = 5-af(x,t)a(x,<), (5) 

S = S/aQ, where operators a^/a are magnon creation/annihilation operators 
satisfying the (bosonic) commutation relation; [a{x,t), {x.' ,t)] = S{x — x'). 
Up to the 0{S) terms, localized spins reduce to a free boson system. 

Consequently in the quadratic dispersion (i.e. long wavelength) approxima- 
tion, the localized spin with the applied magnetic field along the quantization 
axis (z-axis) is described by the Hamiltonian Hmag! 

Hmag= / dxat(x,i)(-^ + s)a(x,0, (6) 

Jxg(intGrfacG) ^ ^™ ' 
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and the Hamiltonian, TLcxi— + ^cx)j '^^^ be rewritten as 

= -Jsj dxc^K,t)a'cijc,t), (7) 

J x£ (interface) 
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G (interface) 



45 



1 - 



a' (x, t)a(x, t) 



AS 



a(x,t)c^(x,t)CT"c(x,t) L (8) 



The variable m represents the effective mass of a niagnon. We have denoted a 
constant appUed magnetic field along the quantization axis as B, which includes 
g-factor and Bohr magneton. In this paper, we take h = 1 for convenience. 

We have adopted the continuous limit and hence, as shown in Fig. [1] (b), 
the interface is defined as an effective area where the Fermi gas (conduction 
electrons) and the Bose gas (magnons) coexist to interact. The dynamics is 
described by the Hamiltonian H'^^ and the width of the interface is supposed to 
be of the order of the lattice constant. [HI Eq. ([5]) shows that localized spins at 
the interface lose spin angular momentum by emitting magnons and conduction 
electrons flip from down to up by absorbing the spin angular momentum (see 
Fig. H]), and vice versa. This Hamiltonian H'^^, which describes the interchange 
of spin angular momentum between localized spins and conduction electrons, is 
essential to spin pumping mediated by magnons. 

Therefore we clarify the contribution of magnons accompanying this ex- 
change interaction to spin pumping, under a time-dependent transverse mag- 
netic field. This is the main purpose of this paper. We also discuss how to 
enhance the spin pumping effect. 

This paper is structured as follows. First, through the Heisenberg equation 
of motion, the spin relaxation torque, which breaks the spin conservation law 
for conduction electrons, is defined in ^ Second, we evaluate it through the 
Schwingcr-Kcldysh formalism in ^ Third, we discuss how to enhance the spin 
pumping effect mediated by magnons through quantum fluctuations in ^ Last, 
we discuss the distinction between our quantum spin pumping theory and the 
one proposed by Tserkovnyak et al. [2] in ^ 

Let us remark that in this paper, we use the term, quantum fluctuations, to 
indicate applying time-dependent transverse magnetic fields. 

2 Spin Relaxation Torque 
2.1 Theoretical model 

Wc apply a time-dcpcndcnt transverse magnetic field, which acts as a quantum 
fiuctuation; 

r{t) EE Tocosint). (9) 
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Then the total HamHtonian of the system (mterface), is given as 



n = H„,ag + H;, + Hci + Ki+Kna.> where 



(10) 



f r B 
Hci - / dxct(x,i) (J5'+-)a" c(x,t) 



(11) 



J x£ (interface) 



Ki = ^ / dxct(x,t)(a++a-)c(x,0. 



(12) 



(interface) 




(13) 



The variable md denotes the effective mass of a conduction electron. 

Here let us mention that in the last section, we have rewritten localized spin 
degrees of freedom into magnon ones via the Holstein-Primakoff transformation. 
There, of course, off-diagonal terms (magnon-magnon interactions) may emerge 
as well as eq. ©, but they cannot satisfy the resonance condition discussed at 
^5.11 in detail. Consequently, their contribution to spin pumping is extremely 
smaller than V^^^ (eq. (|64| ) and negligible. Therefore we omit such interactions 
and concentrate on the effect of V^^agj which represents the coupling with the 
time-dependent transverse magnetic field. 

2.2 Definition 

The spin relaxation torque (SRT),[7] , is defined as the term which breaks 
the spin conservation law for conduction electrons; [HI IS] 



where the dot denotes the time derivative, js is the spin current density, and 
pI represents the z-component of the spin density. We here have defined the 
spin density of the system as the expectation value (estimated for the total 
Hamiltonian, "H); 



In this paper, we focus on the z-component of the SRT. 

Through the Heiscnberg equation of motion, the z-component of the SRT is 
defined as 



(14) 



(15) 




(16) 
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This term arises from 'H^^ and V^^, which consist of electron spin- flip operators; 

rs^ = [Ps,^cx + Ki]A- (17) 

Eq. shows that the SRT operates the coherent magnon state. [ID] 

According to Ralph et al., the net flux of nonequilibrium spin current 
(i.e. the net spin current) pumped through the surface of the interface, J • 
rfSinterface, Can be Computed simply by integrating the SRT over the volume of 
the interface; 

[ d^Z'= [ dSintcrfacc (18) 

'^xG{intcrfacc) 

In addition in our case, conduction electrons cannot enter the ferromagnet, 
which is an insulator. [3] Therefore as shown in Fig. [T](b), the z-component of 
the net spin current pumped into non-magnetic metal can be expressed as 

d^V- (19) 

G (interface) 

Thus from now on, we focus on 7^^ and qualitatively clarify the behavior 
of the spin pumping effect mediated by magnons at room temperature. Let 
us mention that the above relation between the SRT and the pumped net spin 
current can be understood via the spin continuity equation, eq. (|14p . and we 
discuss at i J4.1l 



2.3 Magnon continuity equation 

Here let us emphasize that the spin conservation law for localized spins (i.e. 
magnons) is also broken. The magnon continuity equation for localized spins, 
[12] which corresponds to the equation of motion for localized spins [T3j and 
describes the dynamics, reads 

Pm + v-jr„ = (20) 

where jm is the magnon current density, and represents the z-component of 
the magnon density. We have defined the magnon density of the system also as 
the expectation value (estimated for the total Hamiltonian, "H); 

pf„^(ata). (21) 
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In addition, we call 7^ the magnon source term [T^] , which breaks the magnon 
conservation law. This term arises also from H'^^ and V^^^^^; 



' m Lrm' ' ^cx ' *^magJ/ 



(22) 



aT(x, t)a(x, t) 



AS 



c^{x,t)(T+c{x,t) 



1 - 



a^(x, t)a(x, t) 



AS 



a(x, i)c^(x, t)CT c(x, t) 



a^(x, t)a(x, t) 



AS 



a(x, — a^(x, t) 



1 - 



a'''(x, t)a(x, i) 



45' 



Within the same approximation with the SRT (, which is discussed at the 
next section in detail), this magnon source term, in fact, satisfies the relation; 



(24) 



Then, the z-component of the spin continuity equation for the total system (i.e. 
conduction electrons and magnons) becomes 

Ptotal+V-jLtal -0, (25) 

where the density of the total spin angular momentum, p\^^^\i is defined as 

Ptotai = - Pm, (26) 

and consequently the z-component of the total spin current density, jtot^i' ^"6- 
comes 



Jtotal — Js Jm (27) 

(note that, S^ = S — a^a, via the Holstein-PrimakofF transformation). The spin 
continuity equation for the whole system, eq. (|25p . means that though each spin 
conservation law for electrons and magnons is broken (see eqs. and ((20| ). 
the total spin angular momentum is, of course, conserved. 



3 Schwinger-Keldysh Formalism 

The interface is, in general, a weak coupling regime: [H] the exchange interaction, 
J, is supposed to be smaller than the Fermi energy and the exchange interaction 
among ferromagnets. Thus H[,^ can be treated as a perturbative term. In 
addition, we apply weak transverse magnetic fields. Then we can treat T-L'^^, 
V^y, and V,]iag perturbative terms to evaluate the SRT, eq. pJ|) . 

Through the standard procedure of the Schwinger-Keldysh (or contour- 
ordered) Green's function (see also Appendix IA.2[) .p^ [TCI [T7 ] the Langreth 
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method (see also Appendix ET1),[I1 [H [201 HI] each term of the SRT can be 

evaluated as follows; the first term of eq. reads (see also Appendix IA.3[) 



iJal\l-a\^,t) 1 
2^2> J (2^)3 



at(x, t)a{x, t) 



AS 



cLu 
2^ 



1 - 

S 



c\x,t)<7+c{yL, t) 



(27r)3 



do;' 



(28) 



-2iOf <^ a 



Ga I y^a 



-^0{J^)+0{T^)+0{JS~^) 



The variable tj'(<'>) is the fcrmionic time-ordered (lesser, greater) Green's func- 
tion, and G^'''' is the bosonic lesser (advanced) one. 

We here have taken the extended time (i.e. the contour variable) defined on 
the Schwinger-Keldysh closed time path, [211 [13 El [HI [2Q] c, on the forward 
path C-i. (see also Fig. A-1 in our manuscript [12] ) : c = c^ +c^. Even when the 
time is located on the backward path c<_ , the result of the calculation does not 
change because each Green's function is not independent; Q'^ — Q'^ = Q'- — ^ 
where Q'^^^^ represents the retarded (advanced) Green's function. This relation 
comes into effect also for the bosonic case (see eq. ([79])). 

Here it would be useful to mention that under the thermal equilibrium condi- 
tion where temperature difference does not exist between ferromagnet and non- 
magnetic metal, the ©(J^F*^) term including no quantum fluctuations cannot 
contribute to spin pumping because of the balance between thermal fluctuations 
in ferromagnet and those in non- magnetic metal. [22] 

The second term of eq. ([T6| reads (see also Appendix IA.3|) 



a^"(x, i)a(x, t) 




a{x,t)c^ {x,t)a c(x, t) 
dk' f duj' 



t,k,w-I2»4„k,aj 



^t,k,w-0^i,k,w 



-2intQr 



0,-0 



Go n 



»t:k,w+r!y4.,k,tj »t:k,w+0»i,k,i 



(29) 



where the variable G"^ is the bosonic retarded Green's function. The third term 
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of eq. ([TC)) reads (see also Appendix IA.3 



dk f duj 
2^ 



2 2 



(27r)3 
xJ(e2'^^* + l)GS,_ 



1 - 

5 



dk' 

(2^ 



2TT ^k'--' 



(30) 



+ (e- 



-2int 



l)Ga 



o,si 



^I,k,w+f2^t.k,w ^i^k,w+f2^^k,w 



+C'(J°) + ©( J2) + 0{T^) + C'(JS'-i). 



We omit the ©(J") terms because they contain no contributions of magnons via 
the exchange interaction and are not relevant to the spin pumping effect; the 
terms are out of the purpose of the present study. 

The last term of eq. reads (see also Appendix IA.3[) 



r(i)ct(x,t)cr-c(x,t) 



(27r)3 J 2^ 



1 



dk' 

5 7 ("2^ 



duj' , 
"2^ 



-G 



k' .uj' 



X (e 



t,k,w-n!^i,k,tj 



7< 



^t,k,w-f2!^4,,k,w 



(31) 



+ (e- 



-2iQ.t 



+0(7°) + 0(7^) + ^^(r^) + 0{JS-^). 
Finally, the SRT can be rearranged as 

J2inVLt 



n=0,±l 



where 



£0,2 f dk 
2^2> J (27r)3 



du 
2^ 



dk' 



SJ (27r)3 



dw' 
27 

> 



(Gqji + Go^_f2)(fJ-(, k,aj-0^i,k,Lj ~ ^t,k,Lj-J2^i.k.w 



(32) 



(33) 
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='i.k.Lj+0»t,k,w-' 



(34) 



dk 

(2^ 



did 

2^ 



1 - 

S 



dk' 

(2^ 



diu' 

27^'^''' 



Ga-n)(^lk,w-i2^^t,k,w ~ ^4^k,w-o^t^k,w) 



G, 



t,k,w-nyi,k,w 



^1%k,Lj-J2^i^k,a;) 



,f2 + Go,o)(^t 



* k,tj+J2^I,k,w ^t!k,tj+J2^i!k,Lj3 



(35) 



Here let us mention that, in real materials, there does exist impurity scat- 
tering. We assume that this is the main cause for the finite lifetime of magnons 
and conduction electrons. Moreover the rate of impurities such as lattice de- 
fects and nonmagnetic impurities is, in general, far larger than that of magnetic 
impurities. Therefore we phenomenologically introduce the lifetime and regard 
it as a constant parameter. Then we adopt Green's functions including the the 
effects of impurities as the lifetime, and calculate the SRT by using them; eqs. 

dSSli-dlSl). 

Though, as the result, it might be better to execute the accompanying vertex 
corrections from viewpoints of theoretical aspects, we have not done for the aim 
now explained; to put it briefly, in order to clarify that pumped spin currents 
are generated purely by quantum fluctuations, we have not executed vertex 
corrections. 

It should be noted that, before our present study, Takcuchi et al.[23] have 
already studied spin pumping, on the basis of Schwingcr-Kcldysh formalism, 
under the same condition with ours except two points; (i) they have treated 
localized spins as not magnons but classical variables, and (ii) they have not 
applied any transverse magnetic fields. On their condition, they have clarified 
that, under the uniform magnetization, spin currents cannot be generated with- 
out vertex corrections (i.e. multiple scatterings of impurities). In other words, 
they have already revealed that spin currents can be generated by the effects of 
multiple scatterings of impurities, i.e. vertex corrections. 

Thus, the main purpose of the present study is to propose an alternative way 
for the generation of spin currents without using vertex corrections, i.e. multiple 
scatterings of impurities; we propose a method for the generation of spin currents 
by using time-dependent transverse magnetic fields, which are under our control 
and act as quantum fluctuations. Therefore we call this method quantum spin 
pumping. In order to clarify that pumped spin currents are induced purely by 
quantum fluctuations, we have not included the effects of multiple scatterings 
of impurities (i.e. we have not executed vertex corrections). 



10 



Each Green's function including the the effects of impurities as the lifetime 
reads as follows (see also Appendix IA.2I) : [TB] 

K-Wk'-*/(2T,„)]-i (36) 

(GL',.')*, (37) 

[c^-c^a,k- j/(2r)]-i (38) 

(^^;,k,J*, (39) 

-/b(^')(G^',..'-GL^.O, (40) 

/F(c^)(t;:,k,.-e;.k,J, (41) 

Q'.M.^+5aM.u (42) 

GIk^+^I.Ku. (43) 

where the variables Tm, t, /g(w'), and /f('^) arc the lifetime of magnons, that 
of conduction electrons, the Bose distribution function, and the Fermi one, 
respectively. The energy dispersion relation reads Wk' = Dk''^ + B and Wo-^k = 
Ffc2 - ( JS* + B/2)cr - where £> = l/(2m), F = l/(2mei), a = +1, -l(=t, \), 
and fjL denotes the chemical potential. 

Wc consider a weak magnetic field regime and omit the 0{\{JS + B)/eYY) 
terms, where ep represents the Fermi energy. Through the Sommerfeld expan- 
sion, the chemical potential is determined as, /^(T) = ep — (7rfcBT)^/(12eF) + 
0{T^). 



Gk'.tj' = 
Q^v = 

^(T,k,t*j 



4 Spin Pumping 

4.1 Breaking of spin conservation law 

The relation between the SRT and the pumped net spin current in ij2.2[ /^g (interface) 
J Js -dSintcrfaco, cau be understood via the spin continuity equation; + V -jf = 

' s ■ 

Through the same procedure and the same approximation with the last 
section, the time derivative of the spin density is estimated (on the condition 
mentioned at §5.ip as, p^/T/ ^ 10~^. Thus in the spin continuity equation 
is negligible [121 ^1^'^ our case. 

Consequently, the spin continuity equation becomes V • jg = 7^^ . Therefore 
by integrating over the volume of the interface, the relation between the SRT 
and the pumped net spin current mediated by magnons is clarified; 

/ dxV-jf. (44) 

'J (interface) 

J Jb ■ (^Sinterface- (45) 

In addition in our case, conduction electrons cannot enter the ferromagnet, 
which is an insulator. Then the net spin current pumped into the non-magnetic 



dx = 

xG(intcrfacc) 
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metal can be calculated by integrating the SRT (eq. (pSj) ) over the interface (see 
Fig. [1} . This is our spin pumping theory mediated by magnons via the exchange 
interaction at the interface between a non-magnetic metal and a ferromagnetic 
insulator. The breaking of the spin conservation law for conduction electrons, 
Ts ^ is essential to our spin pumping theory. 



4.2 Quantum fluctuation 

Our calculation (eqs. ([5^ - ([55|) ) gives, 

r/ 0. (46) 

Thus our formalism (eq. (1451) ) shows that spin currents mediated by magnons 
cannot be pumped without quantum fluctuations; this is the significant differ- 
ence from the theory proposed by Tserkovnyak et al. [5] , which is discussed at ^ 
in detail. That is, quantum fluctuations are essential to spin pumping mediated 
by magnons as well as the exchange interaction between conduction electrons 
and ferromagnets (see eqs. ([5^ - ([55|) ): 

V « ^r^. (47) 

This is the main feature of our quantum spin pumping theory mediated by 
magnons. In other words, we have shown that magnons accompanying the 
exchange interaction cannot contribute to spin pumping without quantum fluc- 
tuations. 



4.3 Resonance 

The SRT (cq. ([32])), , is rewritten as (see also Appendix E| 

V = a / (\/—dfc) r/, where (48) 

Jo V ep 

V ^ Yl V{n)e'^-''\ (49) 

n=0,±l 

2 



A = . (50) 

4(27rVF)3 ^ ' 

The variable, 7^^, is the dimensionless SRT as a function of the wavenumber for 
conduction electrons y^F/epk, the angular frequency of the applied transverse 
magnetic field fi/ep, the magnitude of the exchange interaction J/ep, and so 
on. Fig. [2] (a) represents the 

^^(O), (51) 
which is the time average of 7^^ . Fig. [2] (b) shows 

r//A (52) 
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Figure 2: (Color online) (a) The angular frequency dependence of the SRT; 
TJ^iO). A sharp peak exists on the point where the resonance condition, Q ~ J, 
is satisfied. Each quantity, '7^^(±1), has the same structure with T/iO). (b) 
The time evolution of the SRT, at the resonance point {ft = J); (i) 

0.446 X cos2(f7ot), (ii) 0.112 x cos2(f]oV2), (iii) 0.0281 x cos'^ {not / A) , where 
rio = 1.70 X 10^3 



(see eq. (05])), which describes the time evolution of the SRT. We have set each 
parameter, as a typical case,[T31 [MJ [3J [5S] as follows; ep = 5-6 eV, B/ep = 
1.2 X 10-^ T = 300 K, = 4 eV A^, D ^ 0.3 eV A^, t 1 ps = lO^^ x r^, 
ao = 3 A, S* = 1/2. 

Fig. [2] (a) shows that the SRT has a sharp peak as a result of the resonance 
with the angular frequency, il, of the applied transverse magnetic field. The 
sharp peak exists on the point where the condition, 

n = J, (53) 

is satisfied. This is because according to eq. ([7]), the localized spin acts as an 
effective magnetic field along the quantization axis J, which is far larger than 
the applied magnetic field B; 

n^=^/^ = -jf dxs^(x,t). (54) 

J xG (interface) 
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This fact (i.e. resonance condition) is useful to enhance the spin pumping effect 
because the angular frequency of a transverse magnetic field is under our control. 
In addition Fig. [5] (a) also shows that the stronger the exchange interaction 
becomes, the larger the SRT does. Fig. ^ (b) shows TJ^ /A at the resonance 
point (ri = J). Each period reads 

tt/Q. (55) 



5 Distinction from the Theory Proposed by Tserkovnyak 
et al. 

Last, let us discuss the distinction between our quantum spin pumping theory [25] 
and the one proposed by Tserkovnyak et al.[5] It should be noted that, as has 
been mentioned in their article^ (see §VIII. SUMMARY AND OUTLOOK in 
their article [2]), they have phcnomenologically treated the spin- flip scattering 
processes, which we have regarded as the most important processes for spin 
pumping. Nevertheless, now their spin pumping theory has been widely used 
for interpreting vast experimental results. 1271 1301 131j in particular by exper- 
imentalists. Thus, it would be significant to clarify the difference between our 
quantum spin pumping theory and the one proposed by Tserkovnyak et al. 

5.1 The spin pumping theory proposed by Tserkovnyak 
et al. 

According to the phenomenologicalpl I30| spin pumping theory by Tserkovnyak 
et al. and their notation, [32l [2] the pumped spin current Ig-pump reads 

Is-pump = of V X rii + G«m, (56) 

where the dot denotes the time derivative. We have taken e = 1, and m(x, t) de- 
notes a unit vector along the magnetization direction; they have treated m(x, t) 
as classical variables. The variable is the complex- valued mixing conduc- 
tance that depends on the material; [331 [M] G± = G^'' + iG^J^ . 

5.2 Landau-Lifshitz-Gilbert equation 

The magnetization dynamics of ferromagnets can be described by the Landau- 
Lifshitz-Gilbert (LLG) eq.; 

m = 7Hcff X m + am x m, (57) 

where 7 is the gyro-magnetic ratio and a is the Gilbert damping constant 
that determines the magnetization dissipation rate. Here it should be empha- 
sized that though this Gilbert damping constant, a, was phcnomenologically 
introduced. [35] it can be derived microscopically by considering a whole system 
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including spin relaxation; thus the effect of the exchange coupling to con- 
duction electrons should be considered to have already been included into this 
Gilbert damping term. 

The effective magnetic field is set as 

Heff-(r(t),0,B), (58) 

where Tit) represents a time-dependent transverse magnetic field. The LLG eq., 
eq. (|57)) . becomes 

-Bva^ \ / m^rii^ — rh^m^\ 

7 I Bm^ — Fm^ +0^ m^rii^ — rii^m'^ . (59) 
Tmy / \m^uiy - rii^my / 

5.3 Pumped spin currents based on the theory by Tserkovnyak 
et al. at finite temperature 

Eq. (|59p is substituted into I^.p^mpi ^Q- (|56)) : we include the contribution of 
the Gilbert damping term, which depends on the materials, up to 0{a)] a ~ 
10^3 10^2 for NigiFeig (metal), [H] and a - 10"^ for Y3Fe,50i2 (insulator), [3] 
as examples. Their theory is applicable to both ferromagnetic metals and 
insulators. [57] 

Consequently, the z-componcnt of the pumped spin current reads 
If-pump = ^ |7B[(m^)2 + [xny f] - 7rm^m^ - a7r[m^m2'm^ + (m?')^ + (m^ 

+ G^J^ \a-i{B[{uff + (m«)2] - rm^m"} + 7rm^} + 0{a^). 
[G<l"' + aG^']7S[(m-)2 + (m^)^]. 
Gf ^7B[(m")2 + [uvyf]. 



At finite temperature, the magnetization is thermally activated; m 7^ 0.[24) 
Then the time derivative of the z-component means 

rii^ = -iTviiy + a-i{B[{ijff + {myf]-Tm'ui'-} + 0{a^). (63) 

a-tB[{nff + {myf\. (64) 
^ 0. (65) 

5.4 Distinction 



Eqs. (|61]) . ((62)) and (|65p mean that, within the framework by Tserkovnyak et 
al. with the LLG eq., they may gain spin currents at finite temperature if only 
the magnetic field along the z-axis, -B, is applied; 

Is-pump ^ 0- (66) 
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That is, the spin pumping theory by Tserkovnyak et ah [21 [321 with the LLG 
eq. concludes that spin currents may be pumped at finite temperature without 
time-dependent transverse magnetic fields. 

On the other hand, as discussed in the last section (see eqs. (dHj) and ([TZ)) ). 
our approach based on the Schwinger-Keldysh formahsm gives different resuh; 



That is, our quantum spin pumping theory means that spin currents medi- 
ated by magnons cannot be pumped without quantum fiuctuations (i.e. time- 
dependent transverse magnetic fields, r(t)); quantum fiuctuations are essential 
to spin pumping mediated by magnons as well as the exchange interaction be- 
tween conduction electrons and fcrromagnets. 

This (eqs. (j67p and (j66p ) is the significant distinction between our quantum 
spin pumping theory and the one proposed by Tserkovnyak et al. 

6 Summary and Discussion 

Wc have microscopically studied quantum spin pumping mediated by magnons 
by evaluating the SRT at room temperature. Localized spins of the ferromag- 
netic insulator lose spin angular momentum by emitting a magnon and con- 
duction electrons flip from down to up by absorbing the momentum. Thus our 
formalism contains no phenomenological treatments of spin-flip scattering pro- 
cesses. The SRT breaking the spin conservation law represents the (net) spin 
current mediated by magnons, which is pumped into the adjacent non-magnetic 
metal. 

Through the Schwinger-Keldysh formalism, we have clarified that quantum 
fluctuations (i.e. time-dependent transverse magnetic fields) induce a net spin 
current, which can be enhanced through a resonance structure as a function 
of the angular frequency of the applied transverse field. We conclude that the 
breaking of the spin conservation law and quantum fluctuations are essential 
to quantum spin pumping mediated by magnons accompanying the exchange 
interaction. 

In this paper, we have theoretically introduced and deflned the interface as 
an effective area where the Fermi gas (conduction electrons) and the Bose gas 
(magnons) coexist to interact. Though the behavior of the quantum spin pump- 
ing effect mediated by magnons can be qualitatively captured by calculating the 
SRT, the theoretical estimation for the volume, in particular the width, of the 
interface is essential for the quantitative understanding. Of course the width of 
the interface may be roughly supposed to be of the order of the lattice constant, 
but we consider microscopic derivation (so called proximity effects) is an urgent 
theoretical issue. 

In addition, though we have not executed vertex corrections at the present 
study in order to clarify that pumped spin currents are generated purely by 
quantum fluctuations, calculating the SRT under the effects of vertex corrections 




(67) 



16 



(i.e. multiple scatterings of impurities) as well as quantum fluctuations is a 
significant theoretical issue. Also from viewpoints of theoretical aspects (i.e. 
quantum field theory), it might be better to execute vertex corrections. Thus, 
we would like to tackle this issue in the near future. 

In the present study, we have focused exclusively on the dynamics at the 
interface (J 7^ 0), where pumped spin currents are generated with interchang- 
ing spin angular momentum between conduction electrons and magnons. Then, 
our quantum spin pumping theory microscopically well describes the dynam- 
ics of generating pumped spin currents, in particular, the spin-flip scattering 
processes, which we have regarded as the most important processes for spin 
pumping. This is the strong point of our theory; note that the theory proposed 
by Tserkovnyak et al. has phenomenologically treated spin- flip processes. [2] On 
the other hand, the dynamics of pumped spin currents in the non-magnetic 
metal [J = 0), i.e. how the pumped spin currents flow in the non-magnetic 
metal, is out of the application (purpose) of the present our theory. That is, 
though our formalism has microscopically captured the spin-flip scattering pro- 
cesses at the interface beyond phenomenology, it does not cover the dynamics 
of pumped spin currents; the dynamics under the existence of pumped spin cur- 
rents such as the effect of localized spins and conduction electrons on pumped 
spin currents, which is often indicated as the /3 term, [HI [551 US] is out of the 
application (purpose) of the present study. Therefore, we would like to brush 
up our quantum spin pumping theory to cover the dynamics of pumped spin 
currents in the non-magnetic metal; how the pumped spin currents flow in the 
non-magnetic metal. We recognize that clarifying these issues is significant from 
viewpoints of applications as well as fundamental science. 

Last, we have revealed that magnons accompanying the exchange interac- 
tion cannot contribute to quantum spin pumping without quantum fluctuations. 
Here it should be stressed that in our formalism, the meditation of spin angular 
momentum is restricted to only magnons. Wc will take the effect of phonons 
into account and develop our theory as a more rigorous formalism to reveal the 
microscopic (quantum) dynamics of the magnon splitting.}^ 127] Moreover, we 
are also interested in the contribution of magnons to quantum spin pumping 
under a spatially nonuniform magnetization. 

Acknowledgements 

We would like to thank K. Totsuka for stimulating the study, and G. Tatara for 
reading the manuscript and useful comments. We are also grateful to T. Taka- 
hashi and Y. Korai for fruitful discussion, N. Sago for helpful correspondence 
on numerical calculations, and K. Ando for sending the invaluable presentation 
file, prior to the publication, on the 6th International School and Conference 
on Spintronics and Quantum Information Technology (SpintechS). We also 
would like to thank S. Onoda and T. Oka for crucial comments at the 26th 
Nishinomiya- Yukawa Memorial International Workshop Novel Quantum States 
in Condensed Matter {NQS2011). 



17 



We are supported by the Grant-in- Aid for the Global COE Program "The 
Next Generation of Physics, Spun from Universality and Emergence" from the 
Ministry of Education, Culture, Sports, Science and Technology (MEXT) of 
Japan. 

A Langreth Method 

In this section, we show how to calculate eqs. p8)) -(|3T |) in detail. 

First, we briefly show the Langreth method, [121 [HI [HI [ID] which is useful 
to evaluate the perturbation expansion of the Keldysh ( or contour-ordered) 
Green's function[TCl [T71 [TOl [TB] in subsection [XTl Second, we introduce the 
point of the bosonic Keldysh Green's functions |161 [T7[ [201 [?T] in subsection lA.21 
Last, the detailed calculations of eqs. (P5)) - (|5T|) are represented in subsection 

We omit the label, x, when it is not relevant. 

A.l The Schwinger-Keldysh closed time path; a concrete 
example 

For simplicity here, we consider the perturbative term, V' , 

V' = r{t) J dx[a{x, t) + at(x, t)], (68) 

and evaluate the expectation value of the bosonic annihilate operator, a, as an 
example; 

(a(r)> = (i:,a{T)cM-^ jdT'V{T')]) (69) 
= -I f dx' JjT'T{T')(T,a{x,T)a\x.',T'))+0{T^) (70) 
dx'I. (71) 



Here Tc is the path-ordering operator defined on the Schwinger-Keldysh closed 
time path,[3T] c (see Fig. A-1 in our manuscript [T^ [T7 1 [TOl [T5] ) We express the 
Schwinger-Keldysh closed time path as a sum of the forward path, , and the 
backward path, c^; c = + c<_.IlTl [T71 [THl [H] We take r which denotes the 
contour variable defined on the Schwinger-Keldysh closed time path on forward 
path, C-i.. Even when r is located on backward path, c<_, the result of this 
calculation is invariant because each Green's function, G', G'', G^, G^, is not 
independent ; dZl [H dOl \M they obey, 

G"^ - G^ = G> - G<. (72) 

Note that this relation comes into effect also for the fcrmionic case: jl7l[T6l[20l[T9] 

-g^ = g^ -gi. (73) 
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The integral on the Schwinger-Keldysh closed time path of eq. ([7T|). X, is 
executed by taking an identity[T7l [20] into account 



dr' = dr' + dr' 



(74) 



as (see also eq. ([75j) in Appendix IA.2P 

/oo 
dt'T{t')[G\t,t')-G<{t,t')]. (75) 

By using the relation, G''{t,t') = G\t,t') - G'^{t,t'), we obtain 

I^i dt'r{t')G'{t,t'). (76) 



oo 



A. 2 Bosonic Keldysh Green's function 

In this subsection, we show the point of the bosonic Keldysh Green's functions. 
The bosonic Keldysh Green's function, G(r, r'), is defined as }16[ [T71 

G(r,T'):=-*(T,a(r)at(r')). (77) 

Depending on the points where t and r' are located on the Schwinger-Keldysh 
closed time path (i.e. c = c_^ + c^), the bosonic Keldysh Green's function is 
expressed as[T6l [T71 [20] 



'G<{t,t') = -i{a^t')a{t)) , when t e c^,t' e c^ 

G^^/N^ jG>(t,t')--«(a(i)aU^')), when r e c^, r' e c^. 

^^'^'^ ^G*{t,t'), when r e c^,r' e c^. ^ ^ 
^G^(t,t'), when r e c^,t' e c^. 

It should be noted that each Green's function is not independent: |16[ [TTI [20] 

G'-G'' = G>-G<^G' + G* = G> + G<. (79) 

In addition, these relations [50] [5T] would be useful on calculation; 

G* = G^ + G>, (80) 

G* = G<-G% (81) 

G' = G*-G< (82) 

= G> - G*, (83) 

(a^(i)a(t)) = iG<it,t). (84) 
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By executing the Fourier transformation, the lesser and greater Green's func- 
tions for free bosons become [TFl [201 



G^',,.' = -/b(^')(G^'x-GL,,^0 (85) 

= -2mfB{uj')6{uj' -uji,,), (86) 

G^,,,, = -2m[l + fB{oj'mu' ^iOi,,), (87) 

Gk',u:' G^,^^, + G^,^, (88) 

= -2Tn[l + 2fB{io')]SiLu' ~LUk') (89) 

= 2iImGi,,_„,coth(/3w72). (90) 

The last one, G^, represents the Keldysh Green's function [1^1 HZ] and the 
relation is called the bosonic fluctuation-dissipation theorem. [21] 

Fermionic Keldysh Green's function 

It would be useful to compare with the (spinless) Fermionic Keldysh Green's 
function, g{T,T'), which is defined as [M [III Hi [H [13 

5(r,T') :=-*(Tec(T)ct(r')). (91) 

Depending on the points where r and r' are located on the Schwinger-Keldysh 
closed time path (i.e. c = c_s. -I- c<_), the Fermionic Keldysh Green's function is 
expressed as[Il [HI [101 [H [IH] 

'g<{t,t') = i{c^t')c{t)) , when T e c^,t' G c^. 

g>{t,t') = -i{c{t)c^{t')), whcnrec^,T'ec^. 

^^^^'^=^g\t,n whenrec^,r'ec^. ^^'^ 

^g^{t,t'), when T e c<_, r' e c^. 

Note that with reflecting the statistical properties, the sign of the lesser Green's 
function is opposite from the bosonic case. In addition, they satisfy the relation: [l6l 

[H [201 [m [E] 

g' -g^ = g> -g< ^g' + g'' = g> +g<, (93) 



g^ = g' + g< (94) 

= g^' + g^, (95) 

g' = -g' + g> (96) 

= -g^ + g<. (97) 

By executing the Fourier transformation, the lesser and greater Green's func- 
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tions for free Fermions beconie[Tni [201 ttH H] 



G^,. = fFHiGl^^GLJ (98) 

= 2TrifF{uj)S{uj - Wk), (99) 

g>^ = -2Tn[l- fF{io)]S{io~Lu,,), (100) 

Gt. ^^^,. + ^k>. (101) 

= -2TTi[l - 2fFi0j)]SiLO - LOk) (102) 

= 2iIme;^_^tanh(/3a;/2). (103) 



The last one, Q^, represents the Keldysh Green's function [121 HZ] and the 
relation is called the fermionic fluctuation-dissipation theorem. [211 HI] 



A.3 Detail of the calculation; eqs. (I28l)-(I3B 

In this subsection, wc will briefly show how to calculate eqs. ([^5)) - PT|) . 

By adopting the Wick's theorem. [TTll^ the left-hand side (LHS) of eq. 
reads 

(The LHS of eq.dMD) = J dx' J d^" j dr' J dr"r(r')r(T") (l - (^c at(x^rOa(x 

X (Te a(x",r")at(x,r))(Te c^(x, r)c|(x', r')) (Te ct(x', r')4(x, r)). 

(104) 

By employing the relation; 



dr' / dr" + dr' dr" 

dr' [ dr" + [ dr' [ dr" , (106) 



and the Langreth method [HI [El [20] with eqs. (|78p and ([92]), the right-hand 
side (RHS) of eq. (|104p can be expressed as 

(The RHS of eq.dlOl) = ^^-^ / dx' / dx" / dt' / dt"r{t')T{t")(l - -.G<{t',t') 

2 J J J -oo J -oo ^ '5 

X [G\t'\t) - G>{t\t)][gl{t,t')g\{t',t) - g<{t,t')g>{t\t)]. 

(107) 

We also here have taken r on forward path, c^ . As discussed in the last subsec- 
tion, even when r is located on backward path, c<_ , the result of this calculation 
is invariant. 

Here it should be noted that 

G\t'\t) - G>{t",t) = G^{t",t). (108) 
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Then, eq. (|107p can be rewritten as 



(The RHS of eq.dMD) = / dx' / dx" / dt' / dt"r{t')T{t")(l - iG< {t',t" 



' — OO ■/ — OO 

X G^(t", t) [Glit, t')g\ it', t) - g< (t, t')g> (t', t)] . 

(109) 



By executing Fourier transformation, we obtain the RHS of eq. 

Through the same procedure, remained terms are evaluated as fohows; 

(The LHS of eq.®) = J dx' J dx" J dr' J dr"r(r')r(r") (l - ^^x', r-)a(x-, rQ) ■ 

X (Te a(x,r)at(x",r"))(Te ct(x, r)4(x', r')> (T^ c^(x', r')ct (x, r)) 



'^''"^ fdx' fdx" r dt' r dt"T{t')T{t")(i-iG<{t',t'\ 

J J J — OO J — OO S 



2 

^ J J J —OO J —OO 

X G^{tX)[Q\{t,t')gl{t',t)-g<{t,t')g>{t',t)] 



(110) 

' (Teat(x',r')a(x',r')) 



(The LHS of eq.®) = j dx' J dx" jdr' j^dT"V{T)T{T"){l - 



S 

X (Tc a(x',T')at(x",r"))(Tc c^(x, t)c|(x', r')) (Tc c^{^ ,t')c\{x,t)) 

= f dx.' I dx" r dt' r dt"T{t)T{t")(i - iG<{t',t')^ 

J t/ i/ — OC^ 'J — OO S 

X G\t', t")[gl{t, t')g\{t', t) - g< (t, t')g> {t', t)]. 



(HI) 

Note that we have adopted the relation; G' — G< = G'' = G^ — G'. 

(The LHS of eq.®) = J dx' J d^' jdr' jdr"nr)V[r")(l (Tq «t(x-,rO«(x-,rO) 



X (Tc a(x",r")at(x',r'))(Te ct(x, r)4(x', r')) (T^ c^(x', r')c|(x, r)) 
I dx' I dx" I dt' I dt"T{t)T{t")(l^ iG<{t',t')'^ 



JalS r r poo 

^ J J J —OO J — OO S 



X G^it", t')[g\{t, t')glit', t) - g<{t, t')g> {f, t)]. 

(112) 

Note that we have employed the relation; G< — G* = G^ = G' — G>. 
By using Fourier transformation, we obtain the RHS of eqs. ([29| - pip . 



A The Form of the Spin Relaxation Torque 

Finally, let us write down the explicit form of the SRT. We have taken h ^ 1 
and have introduced the dimcnsionless variables as follows; O = O/ep, except 
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T = fceT'/ep and t = e-pT/{h). The dimensionless wavenumber for conduction 
electrons has been rewritten; \J~Fjepk = y. 
The SRT (i.e. eq. (08])), 7;^ reads 

V - (I) + (11), (113) 
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O(j2)+O((ro)") 
J du! dyy'^ 



ct(x,i)cr+c(x,t) ) + 



45 

i_ r dk' r du' 

^ S J (27r)3 i 2^ 



/ rft)ct(x,t)g+c(x,^) ' 
\ 4^ , 



./lint 



-2if2t 



i)g;,._ 
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